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C*^ ' Abstract 

(N ■ 

In this letter we give a systematic derivation and justification of the semiclassical model 
, for the slow degrees of freedom in adiabatic slow- fast systems first found by Littlejohn and 

i-^ ' Flynn [5] . The classical Hamiltonian obtains a correction due to the variation of the adiabatic 

subspaces and the symplectic form is modified by the curvature of the Berry connection. We 
show that this classical system can be used to approximate quantum mechanical expectations 
and the time-evolution of operators also in sub-leading order in the combined adiabatic and 
^ I semiclassical limit. In solid state physics the corresponding semiclassical description of Bloch 

rri- electrons has led to substantial progress during the recent years, see [I]. Here, as an illustration, 

we show how to compute the Piezo-current arising from a slow deformation of a crystal in the 
presence of a constant magnetic field. 
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^ . Consider a quantum system with a Hamiltonian H = H{x, — ieV^:) given by the Weyl quantiza- 

^O I tionof an operator valued symbol i?(g,p) acting on a Hilbert spaced = L^(R")(g)Hf = i^(K",Hf). 

Systems composed of slow degrees of freedom with configuration space M" and fast degrees of free- 
dom with state space Hi are of this form. The small dimensionless parameter e <C 1 controls the 
separation of time scales and the limit e ^> corresponds to an adiabatic limit, in which the slow 
and fast degrees of freedom decouple. At the same time e — >■ is the semiclassical limit for the slow 
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p^ , degrees of freedom. Concrete realizations of this setting are the Born-Oppenheimer approximation 

[2] , the semiclassical limit of particles with spin [3 , Bloch electrons in weak fields [TT] and many 
others, see also [1] and references therein. 

In this letter we show that with each isolated eigenvalue e{q,p) of the symbol H{q,p) there is 

^\f , associated a classical system with an e-dependent Hamilton function h{q,p) = e{q,p)+eM{q,p) and 

a modified symplectic form uj{q,p) = ujQ+efl{q,p). Here M and fl depend also on the corresponding 
spectral projection ttq {q, p) oiH{q,p). The correction to the energy M results from a super-adiabatic 
approximation: the true state of the fast degrees of freedom is in the range of a slight modification 
TT of the adiabatic projector ttq = tto{x, — ieVa:). The correction O to the symplectic form is given 
by the curvature of the Berry connection and takes into account the geometry of the eigenspace 
bundle defined by TTQ(q,p). We show that with the help of this classical system one can approximate 
expectations for "slow" observables a ® l^j and also their time-evolution in the Heisenberg picture 
with errors of order ©(e^). 

The classical system described above appeared first in the seminal work of Littlejohn and 
Flynn |5j in the context of WKB approximations. But they neither claimed nor proved the state- 
ments of the present paper. Independently, Niu and coworkers [1] applied the classical model with 



enormous success in the context of the semiclassical description of Bloch electrons. Here the classi- 
cal description provides simple and straightforward derivations of formulas that are hard to justify 
by other means. 

Given the abundance of slow-fast systems in physics, a complete understanding of their semi- 
classical limit is desirable. The main novelty presented in this letter are general and systematic 
proofs showing that the classical model indeed approximates the quantum mechanical expectations 
of time-dependent "slow" observables. Our approach differs from earlier ones in two ways. Based 
on [B] , it is intrinsically gauge invariant as it does not use a local choice of eigenfunctions for the 
fast system. And it directly applies to arbitrary states, not only to semiclassical wave packets. For 
a mathematically rigorous formulation of our results we refer to [^ . 

The structure of this letter is as follows. We first recall a few basic facts about the Weyl calculus 
and the adiabatic approximation. Then we construct the classical Hamiltonian system and prove 
the various semiclassical approximations. Finally we show how to incorporate also time-dependent 
Hamiltonians and compute, as an illustration, the Piezo-current in the presence of a magnetic field. 
We conclude with some remarks on the literature. 

1 Preliminaries and definitions 

1.1 Weyl calculus 

With a function A : M^" — > C{T-Lf) on classical phase space taking values in the linear operators on 
T-Li one associates the Weyl operator 

{M){x) - j^^ I ApAyA {^{x + y),p) e'^^^-yy^i^iy) 

acting on ip £ H = L^(R",Hf). The composition of operators AB = C induces a composition of 
symbols denoted by C = Ajj^B and called the Moyal product. The asymptotic expansion of A^B 
starts with 

A#B X AoBo + e{AiBo + AqBi - ^{^0,^0}) + 0{e^), 

where {•,•} denotes the Poisson bracket {Ao,Bo} = J2^-^i {dpjAodq.Bo ~ dq-Aodp-BoY Since 
Aq and Bq are operator valued functions, they do not commute in general and neither do their 
derivatives. Hence, in general, {A, A} ^ 0, but, if the derivatives of A are trace-class, 

tr({A,A}) = 0, (1) 

because of the cyclicity of the trace. For later reference we state also the subprincipal symbol for 
triple products 

(A#B#C)i = AiBoC^o + AoBiCo + AoSoCi (2) 

-|(Ao{Bo,Co} + {Ao,So}Co + {Ao|Bo|Co}). 

Here and in the following we use the shorthand 

{Ao|So|Co} := dpA^ ■ Bo 9,Co - 9,^o • Sq 5pCo . 



If a symbol A — a ■ id is a. scalar multiple of the identity, then A and all its derivatives commute 
with any _B. As a consequence one can show that in this case 

Aa#Bo - Ba^Aa x -ie{Ao, Ba} + ©(e^) . (3) 

The fact that the remainder term in ^ is of order e^ and not only e^ is at the basis of higher 
order semiclassical approximations. Note that ttq ~ ttq implies that any partial derivative Ojitq is 
off-diagonal with respect to ttq, 

djTTo = 7ro(5j7ro)7r^ + Tr(j-(9j7ro)7ro , (4) 

with tTq = 1 — ttq. Thus for scalar symbols a — al-^f 

{a, ttq} = TToja, tto}ttq + TTo{a, ttoJtto. (5) 

Finally one can express the trace of a product of Weyl operators by a classical phase space integral, 

Tt{AB) = — i— / dqdp lT{A{q,p)B{q,p)) . (6) 

Here Tr denotes the trace on % and tr the trace on %{. 

1.2 Adiabatic approximation 

Let eo(<z,p) be a non-degenerate eigenvalue band of the principal symbol Hf){q,p) of the Hamiltonian 
H = Hq + eHi and 7ro((j,p) the corresponding family of rank one spectral projections. Then there 
exists an associated subspace nH of Ti that is adiabatically invariant: 

Claim: There exists a projection tt with symbol n ~ ttq + etti + 0{e'^) such that [tt, H] ~ 0{e°°). 
As a consequence, the range of tt is almost invariant, [7r,e^'^~] = 0(6°° |i|), and the adiabatic 
approximation e~''^~7r = e~"^^'^~7r -|- C'(e°°|t|) holds. Moreover 

TToTTlTTo = ^ TToItTo, TToItTo . (7) 

Proof. The construction of tt was first done in [6] and is by now standard, c.f. [4]. Equation 
([7]) simply follows from the fact that tt is a projector in the Moyal algebra, = tt^^tt — tt = 
e(7ri7ro + ttotti - ^{ttq, tto} - tti) -1- 0{e'^). 

1.3 The classical Hamiltonian 

It therefore suffices to study the restriction ttHtt of H to the adiabatic subspace itH in order to 
understand the dynamics within ttH. We will show that its semiclassical limit is governed by the 
associated scalar Hamilton function 

h := e + eM, 

with e = eo + etr{Hi7TQ) and Af := ^ tr ({7ro|i/o|7ro}). 

Claim: It holds that 7r#/i#7r - 7r#7f#7r = 0{e^) and thus 7r/i^-7i-if7i- = 0{e^). 

Proof: In the expansion of tt^/i^tt — -k^Hjj^-k the principal symbol vanishes and the subprincipal 



symbol is Mttq + ^{Tro{HQ ~ cq, ttq} + {ttq, iJo — eo}7ro + {7ro|Ho — eo|7ro}). To see that it also vanishes, 
note that (P and ^ imply 

Mttq = ^trdTToli^okoDTTo = 5tr({7ro|i/o - eolTToDvro 

= 5tr(7ro{7ro|i/o -eo|7ro}7ro)7ro 

= 57ro{7ro|i?o -eo|7ro}7ro = ^{Troli/o - eo|7ro} . 

Moreover = dj[{Ho — eo)7ro)7ro = dj{Ho — eo)7ro + {Hq — eo)djnoTro implies that — {7ro,-ffo — 
eol'^'o = {■"'ol^o — eo|7ro}7ro = {ttqIHq — eo|7ro}, which proves the claim and shows also M = 
^tr {{tto\Ho - CoItto}) = -^tr (ttoIttq, iJo - eo}). 

1.4 The symplectic form and its Liouville measure 

In order to obtain semiclassical approximations up to errors of order e^, one needs to take into 
account that the restriction to the range of tt also induces a modified symplectic form uig on ffi.^" 
given by 

^ ^ _ / E„ \ / f]«« n-jp \ 

where the components of fl in the canonical basis are 

^Q/3 := -itr-Hf (7ro[az„7ro,52^7ro]) , 

with z = {q,p) and a, /3 — 1, ... , 2n. By definition f2 is skew-symmetric and one readily checks that 
it defines a closed 2-form. Actually ft is the curvature 2-form of the Berry connection. For e small 
enough We is thus a symplectic form. The Liouville measure Aj associated with the symplectic form 
We has the expansion 

n 

^^ = (l + lE(^L''-^'l))d9'A---Adp" + 0(e2) 
= (1 + ietr (ttoItto, no})) dq^ A ■ ■ ■ A dp" + 0{e^). 

2 Results: Semiclassical approximations 

2.1 Equilibrium expectations 

Let a : R^" — s- R be integrable and / : R ^> R be smooth, then 

TT{iTfiH)a) = j^ J dK fihiq,p))a{q,p) + 0{e'-n- (8) 

Proof: In the following computation we use that [fr,!!] = 0{e°°) and [tt, /i] = 0{e) imply also 
[ttJ{H)] = 0(e°°) and [nj{h)] = ©(e), and, 7r/(H)7r = /(ttHtt) + 0{e°°) and TTf{h)n = 



7Tf{TThTr)TT + 0{e^). Modulo ©(e^-") we find 

TT(TTf{H)(i] = TT(Trf{H)Tra] =Tr(TTf{TTHTT)Tra 

= Tr I 7r/(7r/i7r) 7ra] := Tr I 7r/(ft,) tt a 

Next note that for scalar symbols the functional calculus for pseudo-differential operators implies 
that f{h) - f{h) = 0{s^). Hence with © we have up to ^(e^-") 



Tr 



{iTfiH)a 



dqdp 



(27re)" 



tr{f{h{q,p)){7T#a#7T){q,p)) 



dqdp 



(27re) 



-/(/i(g,p))tr((7r#a#7r)(9,p)). 



Using ([2]) and the fact that a is scalar, we have that the expansion of A := n^a^TT starts with 
Aq — TToao and 

Ai = TToaiTTo + ao7ri7ro + aoTToTTi 

-^TTojaojTTo} - 5{7ro,ao}7ro - ^aoJTTojTTo} . 

Taking the trace we get up to 0{e^) 

tr(7r#a#7r) = ao + e (tr (ttoAitto) + tr (Tr^AiTTfj-)) 
— flo + eai + e2aotr (ttottitto) 
= a(l + ietr(7ro{7ro,7ro})) , 

where we used ^, tr-^f ({ttojTTo}) — and d?]). 

2.2 Egorov theorem 

Let (/)* be the Hamiltonian flow of h with respect to the symplectic form Wg. Then the Heisenberg 
observable A{t) := e'^? ae"'^? can be approximated by transporting the symbol a of A{0) along 
the classical flow 6i, 



n [A{t) 



n = 0{e^) . 



(9) 



Proof: With a{t) := a o 0* we need to show that ^7ra(i)7r = -[7rHn,Tra{t)n] + 0{e^). Let /i2 := 
£~^(7r#iJ:^7r — TT^h^Tr), then 7r:^/i2#7r = /i2 + 0{e°°) and hence its principal symbol satisfles 
(/i2)o = TTo (/i2)o ttq. Thus 



Tt{H — h)TT, TTa{t)TT 



le 



h2,7r#a{t)i^7r + 0{e^) 



is of order e"^ and 



TtHtT, TTa{t)TT 



7ThTT,TTa{t)TT + 0{e'^) 



h,a{t) n - -TT [hn^ait) - a{t)Tr^h) TT + 0{e^) 



h, a{t) 



[h,Tr], [a{t),Tr] tt + 0{e'^) . 



Since h and a{t) are scalar, ^ implies 

i[/i,^)] = Op'^{{h,a{t)}) + 0{e') 

i[^),*] = OpW({a(i),^o}) + 0(£'), 

where Op^(-) := (•), and thus, again modulo ©(e^), 

-[ttHtt, 7ra(i)7r] = 

= TT Op* (j/^, a(i)} - ie[{h, no}, {a{t),7To}]) n 

= TT OpW {{h, a{t)} - ietr {Tro[{h, ttq}, {a(t), ttq}])) ^ • 

In order to compute dta{t) recall that the Hamiltonian vector-field Xh with respect to tu^ is X^ 
~'{uje)°'^dph and thus by definition of a{t) and fi we have 



da{t) 
dt 



dqa 
dp a 



-e flPP En + e^Pi 



dqh 
dph 



= {h,a{t)} -ietr^j (7ro[{/i, ttq}, {a(t), ttq}]) 



2.3 Transport of Wigner functions 

For ^jj G ttT-L wc define the band Wigner function as 

wt, := (1 - i£tr„, (ttoItto, ^o})) tr W^ , 
where W^ is the standard Wigner function. Then wto — wf^ o c/)^* + 0{e'^) in the sense that 



with ^(i) := e ^^ ci/; for all scalar bounded symbols a. 

Proof: Using the invariance of A^ under the Hamiltonian fiow (/)* (Liouville's Theorem), one directly 
computes 



dAa<°C*)«= / dAe<(ao0*) 



= (^, i(i)^)« + 0(e2) = (V.(t), «^W)«+0(£'). 



2.4 Time-dependent Hamiltonians 

It is straightforward to generalize the above statements to the case of a Hamiltonian H{t,q,p) 
depending explicitly also on time. To this end one just adds the canonical pair (i, E) and applies 
the previous results to the symbol K(t, q, E,p) = E+H{t, q,p). Its spectral projections 7ro(i, q,p) are 



independent of E and the classical Hamilton function is k{t, q,E,p) = E + h{t, q, p) with symplectic 
forma; — a;o+er2(t, (7,^), where h and Vi are computed from the instantaneous Hamiltonian i?(i, q,p) 
as before. The equations of motion are now 

''\fl„h-FVlPP8,.h + FnP* 

(10) 



q = 


= (l + el7P«)5p/i-er!PP9,/i + er!P* 


V = 


= (-l + £f^«P)<9,/i-ef}««ap/i-£f}«*. 



On the side of quantum mechanics we have (e '^ ^ V')(^) — U'^{t,to)ip{tQ), where U'^{t,to) is the 
unitary propagator generated by H{t). The statement of the Egorov theorem becomes 



nto)iA{t) - a o 4''°)Hto) = 0{e^) , (11) 

where a o t/)*-*" (g, p) = a{Q{t\to,q,p),P{t\to,q,p)) and Q{t\to,q,p),P{t\to,q,p) is the solution to ([TU| 
with initial data (3(to|io, 9,p) = q and P(io|io, 9,^) = P- 

2.5 The Piezo current in a magnetic field 

As an illustration consider nonintcracting electrons in a slowly deformed periodic crystal subject to 
a constant external magnetic field Bq + eB, where Bq has rational flux per unit cell. The crystal 
is modeled by a potential Vr(i) periodic with respect to some lattice F. After a magnetic Bloch- 
Floquet transformation the Hamiltonian is the Weyl quantization q H> ieVp, p H> p, of the operator 
valued symbol 

Ho{t,q,p) = ^{p-iVy + lBoy+^Bqf + Vr{y,t) 

pointwise acting on L'^{My), where M is the fundamental domain of a magnetic super-lattice F. It 
is convenient to write Bij ~ eijk{B)k for the following computations. Up to unitary equivalence 
HQ{t, k) is periodic in k = p-|- hBq with respect to the dual lattice F* whose fundamental domain is 
the magnetic Brillouin zone T*. To each isolated magnetic Bloch band eo(i, k) of Hq with spectral 
projection TTo{t, k) we associate the classical system h{t, k) = eo{t, k) +eM{t, k). For il{K) one finds 
nP* = -itr(7ro[V^o,ft7ro]), Q^.f - -itr(7ro[9p.^o, 9p,7ro]), QP'^ = ^^ppB, n^i = -{B^ppB. The 
classical equations (|10l) thus read 

q = {l + enPP{t,K)B)d^h{t,K) + enP'{t,K), k = Bq. 

Note that all classical objects are F*-periodic and hence the classical phase space is really R'' x T*. 
We can now approximate the current operator J(i) := ■j^-^(i) on the range of tt (to) according to the 
Egorov theorem (fTT|) and find that the current density when starting in the state p(to) — f{H{to)) 

is 



lim ^ReTr(7:{to)f{H{to))J{t)xAiq)) 
1 f dAe(t) 



A^R3 TA 



e Jf. {2TTf 



{f{hm))o4,Y%n)q{t,K) + 0{e). 



With dAe(i) = (1 + eVL{t) ■ B)An and for a fully occupied band f{h{to,K)) = 1 the leading term 
vanishes as can be seen by partial integration. The contribution to the Piezo current density from 



this band is therefore 



Jit) = l^^j^[Bn{t,K)-Veo{t,K) + np\t,K)) 
■ np\t,K,). 



T' (27r)3 

T. (2Tp 

The last equahty follows again by partial integration and the fact that fl^P is a closed 2-forni 
and thus the divergence of il with (51).; = ^^ijk^^jk vanishes. From this expression it is now 
straightforward to derive the King-Smith and Vanderbildt formula ^ for the orbital polarization, 
cf . [TJ 12] . This shows that the orbital polarization is a geometric quantity that does not change 
under variations of the magnetic field as long as the Fermi energy lies in a gap between (magnetic) 
Bloch bands, see also [TU] for systems with disorder. Note, however, that the result does not follow 
from the modified equations of motion ([TOl) alone, but requires their correct application including 
the modified Liouville measure. 

2.6 Concluding remarks 

As mentioned before, the literature on adiabatic slow- fast systems is vast. Several groups [5|[l][TT| 
[9l 11] arrived independently and with different methods at equations more or less similar to (fTO|) . 
However, their precise connection to quantum mechanical expressions was not established before 
for ([H, and only shown in a special case [IT] for (|9]). 

The most striking applications of the modified semiclassical model are due to Niu et al. p. 
They establish ([TOl) as the equations of motion for the center of a Bloch wave packet. While it 
is natural to conclude from this also the formulas ([U and ([TT]) . they give, to our knowledge, no 
systematic derivation. In particular, they arrive at the correct Liouville measure A^ by looking for 
the invariant measure of (jlO[) and then postulate ([5]). Moreover, in the case of nonzero magnetic 
field Bq the magnetic Bloch bundle defined by 7ro(K;) over the torus T* is not trivializable, which 
might present an obstruction to patching statements about localized wave packets together. Indeed, 
in [11] a rigorous derivation of ([9]) was given for Bloch electrons with Bq = 0, which relies heavily, 
as also [2l[3l|4j[5] do, on the possibility to chose a global non- vanishing section of the Bloch bundle. 
The difficulties with generalizing the approach of [11^ to magnetic Bloch bands led us to the new 
approach presented here. 
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